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^ ' Numerical solutions of Einstein's and scalar-field equations are found for a 

I ■ 

Q^' global defect in a higher-dimensional spacetime. The defect has a (3 + 1)- 



dimensional core and a "hedgehog" scalar-field configuration in n = 3 extra 
dimensions. For sufficiently low symmetry-breaking scales r], the solutions 



X 

I are characterized by a flat worldsheet geometry and a constant solid deflcit 

angle in the extra dimensions, in agreement with previous work. For rj above 
the higher-dimensional Planck scale, we find that static-defect solutions are 
singular. The singularity can be removed if the requirement of staticity is 
relaxed and defect cores are allowed to inflate. We obtain an analytic solution 
for the metric of such inflating defects at large distances from the core. The 
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three extra dimensions of the nonsingular solutions have a "cigar" geometry. 
Although our numerical solutions were obtained for defects of codimension 
n = 3, we argue that the conclusions are likely to apply to all n > 3. 
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I. INTRODUCTION 



According to the "brane- world" picture, our universe is a (3 + 1) -dimensional brane 
floating in a higher-dimensional bulk spacetime [1-3]. The branes can be thought of as 
fundamental D-branes, or they can arise as topological defects in a higher-dimensional field 
theory. It may well be that £)-branes could also eventually find some dual field-theory 
description. In the simplest, codimension-one models, the brane can be pictured as a domain 
wall propagating in a 5D spacetime [1,4]. Higher codimensions with both gauge and global 
defects have also been considered. In particular, models have been discussed where the 
field configuration in the directions orthogonal to the brane is that of a cosmic string [5], 
a monopole [6,7], or a texture [8]. In these models, the bulk curvature produced by the 
defects plays a major role in localizing gravity on the brane and in solving the mass- hierarchy 
problem. 

The solutions of Einstein's equation describing higher dimensional defects with a global 
charge found in Ref. [6] are very simple. In the absence of a bulk cosmological constant, and 
for codimension n > 2, the spacetime is flat along the brane worldsheet and is characterized 
by a constant solid angle deficit in the transverse directions. The angle deficit grows as the 
symmetry breaking scale rj of the defect is increased, and eventually consumes the entire 
solid angle at some critical value rjc, which is comparable to the higher-dimensional Planck 
scale. For rj = rjc, the transverse metric becomes that of a cylinder. 

The solution of Ref. [6] were obtained in the asymptotic region outside the defect core. 
In the core region, analytic solutions cannot be found, and one of the purposes of the present 
paper is to check numerically that the asymptotic solutions of Ref. [6] can be matched to 
appropriate core solutions. 

Another important question is what happens if rj is greater than rjc- One of the mo- 
tivations for looking into this regime is the interesting work by Dvali, Gabadadze, and 
Shifman [9], suggesting that brane- world models of codimension n > 2, with a very low 
value of the bulk Planck scale, can provide a solution to the cosmological constant problem. 
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It was conjectured in Ref. [9] that defect solutions in such models may exhibit de Sitter 
inflation of the worldsheet, with the expansion rate inversely proportional to the brane ten- 
sion. The very low-accelerated expansion rate observed today could then be obtained with 
a very large-brane tension. 

The main purpose of the present paper is to solve Einstein's and scalar-fleld equations 
for super-critical global defects. Specifically, we consider a codimension-3 global monopole. 
In the following Section, we introduce the model and review the asymptotic solutions of 
Ref. [6]. The numerical solutions assuming a static (noninflating) worldsheet are presented 
in Sec. Ill, for both sub- and super-critical regimes. We find that super-critical solutions 
have a curvature singularity at a finite distance from the monopole core. Inflating defect 
solutions are discussed in Sec. IV, where we show that the expansion rate H can be adjusted 
so that the singularity is removed. We interpret the resulting nonsingular geometries as the 
physical defect solutions. Our conclusions are summarized in Sec. V. 



II. FIELD EQUATIONS AND ASYMPTOTIC SOLUTIONS 



The action for our model is 



5 = X\f^^ 



(1) 



where = 1/M2+", M is the 7D Planck mass, and rj is the symmetry-breaking scale. The 
scalar-triplet fleld 0" is assumed to be in a global monopole "hedgehog" conflguration in 
n — 3 extra transverse dimensions, (f)"' — (f>{r)x"'. 

The general static form of the seven dimensional metric with spherical symmetry in the 
extra three dimensions is 



ds' = B'{r)dsi + dr' + C\r)r'dili , 



(2) 



where ds\ represents the 4D world-volume metric. The corresponding Einstein equations 
are 
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where R^"^^ is the Ricci scalar of the 4D world volume. 
The field equation for the scalar field is 



A0(0^ - 77^) = . 



(6) 



The solution obtained in Ref. [6], assuming a fiat 4D world- volume, i?^^) = 0, and (f) — rj 
in the asymptotic region, is of the form 

ds"^ = rji.^dx^'dx'' + dr"^ + (1 - K^rj^ydnl . (7) 

The solid angle deficit is AQ = Ai^K^rf, and the critical symmetry- breaking scale is 

77c = 1/k. (8) 

It was shown in Ref. [6] that for this value of ?], the Einstein equations have a cylindrical 
solution 



ds^ = r]f,^dx^dx'' + dr^ + Cld^ll , 



(9) 



where Cq is an arbitrary constant. The expectation was that this constant radius of the 
cylinder can be determined by matching to an appropriate-core solution, so that the full 
3D geometry will be that of a "cigar". In the next section, we shall see that (i) this 
expectation for the critical solution is not supported by the numerical analysis, but (ii) cigar- 
type geometries do nonetheless make their appearance in super-critical monopole solutions. 
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III. STATIC SOLUTIONS 



We numerically solved the Einstein and scalar-field equations (3)-(6) with boundary 
conditions B{0) = C(0) = 1, B'{0) = C"(0) = 0, 0(0) = 0, and = at large r. For static 
solutions, we assumed a flat worldsheet, dsl — "q^j^udx^dx^ ^ and R^'^^ — 0. 

The solutions for the scalar field and the metric of a sub-critical monopole {kt] — 0.2) 
are shown in Fig. 1. As expected, (j) rapidly approaches 77, and the metric coefficients B and 
C approach constants at large r. The asymptotic value of the angle-deficit coefficient C is 
in agreement with Eq. (7). As 77 gets closer to k~^, the approach to this asymptotic value 
becomes increasingly slow. 

In the critical case, kt] — 1, C ^Oasr— >oo. However, the radius of the extra 
dimensions Cr is a growing function of r (see Fig. 2). Thus, the asymptotic behavior of the 
critical-monopole solution is not described by the cylindrical metric (9). With (f) ^ r], we 
were able to find a power-law analytic solution of Einstein equations, 

^ ^ ^(2-v^)/12 ^ ^-0.097 ^ (10) 
C;^j.{-5+Vl0)/e^j.-0.S06^ (11) 

which gives a good fit to our numerical solution at large r. 

In the super-critical case, kt] > 1, we find that static solutions develop a curvature 
singularity at a finite distance from the defect core. When 77 is only slightly above the 
critical value, this distance can be very large. The solution initially behaves as the "cigar" 
solution (9), but eventually the radius of extra dimensions blows up, C ^ 00, and the metric 
coefficient B plummets to zero at the same finite value of r — rs- This behavior is illustrated 
in Fig. 3 for nrj — 1.06723. The nature of this singularity is similar to that of a static 4D 
global string [10]. 

As the symmetry-breaking scale r] is increased, the singularity gets closer to the core 
of the monopole. We performed numerical calculations for a number of values of r); the 
resulting dependence rs{ri) is plotted in Fig. 4. For rj ^ 1.2, the singularity gets too close to 
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the center, so that the scalar field does not fully relax to its asymptotic value, (t> — rj^ before 
the singularity forms. 



IV. INFLATING SOLUTIONS 



It is well known that the singularity of some static-defect solutions can be removed by 
relaxing the requirement of staticity and allowing the defect worldsheet to inflate. Examples 
are domain walls [11,12] and global strings [13,14] in (3 + 1) dimensions. It is, therefore, 
reasonable to check if the same strategy works in higher dimensions. 

The 4D metric of an inflating defect worldsheet is the de Sitter metric. 



and the 4D curvature scalar is ^^^^ = 12if^. We solved the equations numerically for 
different values of and indeed, for values of r] above, but not too far above the critical, 
we were able to adjust H so that the singularity is removed. Two super-critical, nonsingular 
solutions are shown in Fig. 5. 

We see that the extra dimensions in these solutions have a "cigar" geometry, with 
Cr — > const at large r, and that the scalar fleld quickly approaches a constant value, 
which, however, is somewhat below 77. The latter fact can be easily understood from the 
field equation (6) for 0. With 0' = and \f\r\Cr = Cq = const, it gives 

We were able to find an exact analytic solution of our system of equations, which describes 
the asymptotic behavior of these nonsingular solutions: 



dsl = -dt'^ + e^^^dx^ , 



(12) 



= 00 = 



\/5 — K.'^Tj'^ 



(14) 



VXrjB 



H 

= — sin 




(15) 



where k = ^Jh/l2%{K^rf - 1)/Kr]. 

In order to solve the equations numericaUy, we had to impose the boundary condition 
(f) — (po at large r. The actual sequence of steps that led to the solution was that, after 
experimenting with numerical solutions for different values of H with our old boundary 
condition, 0(r oo) = 77, we guessed that the nonsingular solution should have a cigar 
geometry. We then found the analytic solutions (14)- (16) and used (14) to set the boundary 
condition. 

In the asymptotic solution (15), the expansion rate H can be absorbed by rescaling the 
time coordinate t in Eq. (12); then the solution takes the form 

{X7]^)ds^ = k-\sm^ xdsl + dx^) + Cld^l , (17) 

where x = \/Xrjkr and ds^, is the 4D de Sitter metric with H = 1. The above metric is 
essentially the same as the one found in Ref. [6] as a solution ior (f) — rj and a positive bulk 
cosmological constant. In our model, the role of the cosmological constant is played by the 
scalar-field potential, V{4>q) > 0. 

The solution (17) has an apparent singularity at x = tt, but as it was noted in Ref. [6], 
the first two terms in Eq. (17) describe a 5D de Sitter space. The 5D infiation rate is 
H5 = \/X7]k, and the surface x = ^ corresponds to the de Sitter horizon. This shows that 
the space outside the monopole core infiates not only along the 4D worldsheet, but also in 
one of the transverse directions, while two of the extra dimensions remain compactified in 
a sphere. This behavior is reminiscent of the topological-inflation scenario [15,16], where 
infiation occurs in the cores of topological defects in 4D. There are, however, important 
differences. First, the defect core has a fixed radius and is not expanding in the radial 
direction, and second, two of the extra dimensions remain compactified, also at a fixed 
radius. 

Although the 4D expansion rate H drops out of the asymptotic solution, it is a meaningful 
parameter for the full spacetime. With our boundary conditions at r = 0, it gives the 
infiation rate in the core of the monopole. As we saw in the preceding section, the static 

8 



H — solutions with kt] > 1 are singular. As H is increased, the singularity gets milder, 
and at a certain value of H it disappears completely, resulting in a cigar geometry. We 
interpret the corresponding solution as the physical solution describing the spacetime of a 
superheavy defect. 

Our conjecture is that for each 77, only one particular value of H gives a nonsingular 
solution. Of course, it is impossible to find this precise value numerically, and therefore, 
all super-critical numerical solutions that we considered eventually developed a singularity. 
The nature of the singularity is similar to that of a static solution in Fig. 3, with S ^ and 
C — > 00. By adjusting the value of i7, we were able to shift the onset of the singularity to 
larger and larger radii, so that the solution got closer and closer to the asymptotic solutions 
(14)-(16). In the end, however, C(r) invariably started to diverge, at some point very close 
to the would-be de Sitter horizon, r ~ tt/x/At/A;. We believe that our results are very 
suggestive, but they do not, of course, constitute a proof that all values of H except one 
yield singular solutions. It may be possible to provide a proof using the dynamical-systems 
method employed in Ref. [13]. 

Our numerical results indicate that, for the nonsingular solution, H grows as we increase 
T], and the growth is nearly hnear, at least, for 77 sufficiently close to r]c (see Fig. 6). Above 
some value, k?] > 1.1, it is difficult to get a cigar solution numerically. As r] (H) increases, 
Cr becomes more wiggly in the asymptotic region unless rj and H are tuned very delicately 
to get a cigar solution. This causes numerical difficulties in getting a cigar solution for large 

V- 

From Eq. (14) we see that the asymptotic solution exists only for 1 < ^77 < At the 
end of this interval, Krj — \/5, the scalar field is in the unbroken-symmetry state, (f) — 0. 
Then the source term in the Einstein's equation is pure vacuum, = V{0)6'^, and the 
solution (17) takes the form 

{Xri'^)ds'^ = 8(sin2 xdsl + dx^) + ^d^l . (18) 

The de Sitter horizon radius of the inflating five dimensions in this solution is twice the 
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radius of the remaining compactified two dimensions. Since the field ^ in this solution 
has the same value, = 0, as in the monopole core, the metric (18) describes the entire 
spacetime for kt] = \/5. For higher values of r], we expect that remains equal to zero, and 
the metric is still given by (18). 

V. CONCLUSIONS 

In this paper, we continued the investigation of the gravitational field of higher- 
dimensional global defects that was started in Ref . [6] . We have obtained numerical solutions 
of Einstein's and scalar-field equations for a defect with a (3-|-l)-dimensional core, which has 
a global-monopole-like field configuration in the extra three dimensions. We have verified 
that, for symmetry-breaking scales below the higher-dimensional Planck scale, rj < 1/k, the 
spacetime of the defect worldsheet is flat, and the geometry of extra dimensions is that of a 
generalized cone, in agreement with the asymptotic solution obtained in Ref. [6]. The solid 
angle deficit grows with 77 and consumes the entire solid angle at the critical value rjc — I/k. 

Our main goal in this paper has been to investigate supercritical defects with 77 > rjc- 
First, we have found that static supercritical solutions have a curvature singularity. It has 
been known for some time that singularities of static-defect solutions in 4D can be removed 
by relaxing the staticity requirement and allowing the defect cores to infiate. We have 
shown that the same is true for our 7D global monopoles. We found numerical solutions 
in which the defect core has the geometry of de Sitter space infiating at some rate H. As 
H is increased, for each rj > rjc, there is a value of H for which the spacetime becomes 
nonsingular. The extra dimensions of these nonsingular solutions have a "cigar" geometry. 

At large distances from the core, the nonsingular solutions have a very simple form, and 
we have found an exact analytic solution of the field equations describing this asymptotic 
behavior. In this solution, the five dimensional subspace, formed by the four dimensions 
tangent to the defect worldsheet and the radial dimension in the extra space, has the geom- 
etry of a 5D de Sitter space. Thus, infiation occurs not only along the defect, but also in 
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one extra dimension. (Note, however, that the radius of the defect core is smaller than the 
horizon and remains fixed during the course of inflation.) The remaining two extra dimen- 
sions are compactified as a sphere of a fixed radius. Far away from the core, the scalar field 
approaches a constant value 0o, which, however, is somewhat smaller than the flat-space 
expectation value rj. 

As the symmetry-breaking scale rj is increased, the inflation rate H grows and the asymp- 
totic scalar- held value 0o decreases. For rj > V5k~^, the scalar fleld vanishes in the entire 
spacetime, and the defect solution turns into a vacuum solution of Einstein equations with 
a cosmological term. We interpret this as indicating that, for higher values of r], defect 
solutions with a core of a flxed transverse size do not exist. 

We focused in this paper on a global-monopole-type defect in a 7D spacetime, but since 
the gravitational properties of all global defects of codimension n > 2 are rather universal, 
we expect our conclusions to apply to all such defects, with only minor modiflcations. In 
particular, the analogues of the n — 3 asymptotic solutions (14)- (16) can be easily found for 
n > 3. 

As already noted, the inflation rate H grows as rj is increased. This behavior is opposite 
to what is needed for solving the cosmological constant problem. This indicates that global 
defects cannot be used to implement the Dvali-Gabadadze-Shifman idea [9] . We shall discuss 
the gravity of higher-dimensional gauge defects in a forthcoming paper. 

When this work was in progress, we learned that Ruth Gregory has independently studied 
inflating black-brane solutions. 
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FIG. 1. Scalar and gravitational fields of a sub-critical monopole, kt] = 0.2. The scalar field 
rapidly approaches rj. (The lower part of the graph is not shown.) The metric coefficients and 
approach constants. exhibits the expected solid angle deficit in the asymptotic region. 
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FIG. 2. Metric coefficients B and C for the critical monopole, Krj = 1.0. The asymptotic 
power-law solutions are shown by dashed lines. 
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FIG. 3. Metric coefficients B and C for a super-critical monopole, K-q = 1.06723. The metric 
becomes singular at a finite distance from the center. 
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FIG. 4. The location of the singularity as a function of 77. 
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FIG. 5. Scalar and gravitational fields of the "cigar" solutions with uri = 1.02361 
{H = 0.10543\/A?7; solid line) and kt] = 1.03766 {H = 0.15565\/A7?; dashed line). In the asymptotic 
region, the scalar field and Cr are nearly constant, and B behaves as a sine function. The scalar 
field is scaled to have a better discrimination. 
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